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Abstract. We examine how the initial state (pre-event corona) affects the numerical MHD simulation 
for a coronal mass ejection (CME). Earlier simulations based on a pre-event corona with a homoge- 
neous density and temperature distribution at the lower boundary (i.e., solar surface) have been used 
to analyze the role of streamer properties in determining the characteristics of loop-like transients. 
The present paper extends these studies to show how a broader class of global coronal properties 
leads not only to different types of CMEs, but also modifies the adjacent quiet corona and/or coronal 
holes. 

We consider four pre-event coronal cases: (1) constant boundary conditions and a poly tropic 
gas with 7 = 1.05; (2) non-constant (latitude dependent) boundary conditions and a polytropic gas 
with7 - 1.05; (3) constant boundary conditions with a volumetric energy source and 7 = 1.67; 
(4) non-constant (latitude dependent) boundary conditions with a volumetric energy source and 
7 = 1.67. In all models, the pre-event magnetic fields separate the corona into closed field regions 
(streamers) and open field regions. The CME’s initiation is simulated by introducing at the base of the 
corona, within the streamer region, a standard pressure pulse and velocity change. Boundary values 
are determined using MHD characteristic theory. 

The simulations show how different CMEs, including loop-like transients, clouds and bright rays, 
might occur. There are significant new features in comparison to published results. We conclude that 
the pre-event corona is a crucial factor in dictating CMEs properties. 


1. Introduction 

Coronal mass ejections (CMEs) are recognized as an important component of the 
large-scale evolution of the solar corona and as a key factor in the generation of 
geomagnetic storms. They were first observed with the Orbiting Solar Observatory 
(OSO-7) white-light coronagraph (Tousey, 1973). Later the Skylab coronagraph, 
operating during 1973—1974, recorded 77 events (Munro et al., 1979), and the 
Solwind coronagraph on the P-78 satellite, during 1979-1985, recorded more than 
1200 events of this kind (Sheeley et al., 1980). 

CMEs are observed to occur in a wide variety of sizes and shapes and at various 
latitudes (Munro et al., 1979; Howard et al., 1984). Munro (1977) has classified 
the appearance of mass ejection transients observed during the Skylab period; the 
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most dominant type is the outwardly expanding loop, or loops. Clouds or amor- 
phous blobs constitute the next most common type. The remainder defy specific 
classification. Burkepile and St. Cyr (1993) gave a more detailed description of 
apparent morphologies of CMEs. They emphasize that features often do not fall 
neatly into one category and the shapes of observed features may evolve as they 
move outward through the corona or may be altered due to projection effects if 
they move out of the plane of the sky. Typical characteristics of loop transients 
have been given by Sime, MacQueen, and Hundhausen (1984). 

In CME speed surveys the outward speeds of mass ejection are from less than 
100 km s” 1 to greater than 1200 km s -1 at 1.75 and 6 R (m) , respectively (Rust 
et ai, 1980; Gosling et aL, 1976; Hildner, 1977; Hundhausen, 1977). Recently 
Hundhausen, Burkepile, and St. Cyr (1994) showed that the average CME speed 
variation over solar cycle 21 could be a factor three. This average is low in 1984 
(157 km s _I ), and 1987 (262 km s~ ! ), and high in 1980 (355 km s _l ), 1985 
(458 km s" 1 ), 1986(371 kins" 1 ), and 1989 (410 km s“ j ). 

In order to give a physical interpretation of observed coronal mass ejections, 
numerical magnetohydrodynamic models have been developed, starting in the late 
1970s, first in one and then in two dimensions (Dryer et ai , 1979; Nakagawa and 
Steinolfson, 1976; Nakagawa, Wu, and Han, 1978, 1981; Steinolfson, Suess, and 
Wu, 1982; Steinolfson and Hundhausen, 1988; Wu et ai, 1978, 1982, 1983). In 
these models, the initial conditions, i.e., the pre-event coronae, vary widely. In the 
early models, the pre -event corona was assumed to be a hydrostatic, isothermal 
atmosphere with a magnetic configuration represented either by a dipole potential 
magnetic field (Dryer et ai , 1979; Nakagawa and Steinolfson, 1976; Nakagawa, 
Wu, and Han, 1 978; Wu et ai , 1 978, 1 983) or by a dipolar force-free magnetic field 
in two and half dimensions (Nakagawa, Wu, and Han, 1981 ; Wu etai 1982). Sime, 
MacQueen, and Hundhausen (1984), comparing the observed features of CMEs 
with numerical models, showed that those models are only able to qualitatively 
describe some of the observed characteristics. 

Later models adopted as initial configurations either polytropic steady state 
streamers or quasi-steady-state streamers with an ad hoc energy input (Steinolfson 
and Hundhausen, 1 988, hereafter referred to as SH) and invoked either a thermal or 
a magnetic driving force as the perturbing agent responsible for CMEs’ initiation. 
SH, for example, used three different pre-event atmospheres (hydrostatic state 
with dipole potential magnetic field; polytropic steady-state streamer with dipole 
magnetic field; quasi-steady-state heated streamer with dipole magnetic field) in 
their CME simulations. The structure of these pre-event coronae has been discussed 
by Steinolfson (1988), who suggested that using these modified initial atmospheres 
for CME simulation leads to a better agreement between models and observations. 
In fact, the models of mass ejections which were initiated at the base of the closed- 
field region of a dipole field in a hydrostatic atmosphere failed to simulate the 
most commonly observed features of loop-like CMEs because of the formation 
of a strong fast MHD shock. However, the simulations of a mass ejection in a 
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polytropic corona with a streamer correctly reproduced two features (rarefaction 
between the legs and limited latitudinal leg motion), although the development 
of a fast shock compression front still yielded higher brightness at the loop top 
than in the legs, contrary to what is observed. SH prevented shock formation by 
introducing a volumetric heating term in their energy equation and simulating the 
CME energy source by an increased heating in a localized region within the closed 
field area, thus reproducing successfully most of the loop-like observed features. 
More recently, analogous results have been obtained by Guo ct al (1991), who 
used a similar polytropic streamer without heat source as pre-event atmosphere but 
with an emerging-flux type perturbation. 

In this paper we use four different pre-event coronae, two similar to earlier cases 
and two with more realistically low densities outside the streamers, to further study 
the importance of the initial atmospheres in CME modeling. Section 2 illustrates 
the technique used to construct the pre-event atmospheres. The results of our 
simulations for pre-event coronae are illustrated in Section 3. Their influence on 
the shapes and sizes of coronal mass ejections, as well as the modifications of 
the adjacent coronal hole regions at the time of CMEs events, are discussed in 
Section 4. Concluding remarks are given in Section 5. 

2. Numerical Simulation Method 

The two-dimensional, fully implicit, continuous Eulerian scheme (FICE) in spher- 
ical coordinates (Wang, Hu, and Wu, 1982) has been used to solve the ideal 
magnetohydrodynamic (MHD) equations (Wang et aL, 1993) in the present study 
for both the pre-event coronae and CME simulations. The first version of FICE in 
Cartesian coordinates was given by Hu and Wu (1984). One of the advantages of 
the FICE scheme is to incorporate the characteristic boundary conditions into the 
solutions process for the assurance of self-consistency (Wu and Wang, 1987). 

The computational and physical boundary conditions we adopt for the present 
numerical simulation are the following. 

2. 1 . Computational domain and grids 

In the latitude direction 9 changes from (— A0/2) degrees (near the pole) to (90° + 
A0/2) degrees (near the equator). The spacing in 0-direction is A0 = 4.5° so that 
there are 22 grid points in latitude; the first grid point is at -2.25° and the last 
grid point is at 92.25°. The pole (9 = 0°) and the equator (0 = 90°) are centered 
between grid points. 

The spacing in the r -direction is Ar,; = A0) an dr, = 1 + A0)^ 

The computational domain in the radial direction includes 27 grid points and 
extends from the solar surface to 7 7? . , With this grid selection, the two-dimensional 
grid cells have approximately square shape, which leads to a more evenly weighted 
numerical calculation. 
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2.2. Computational boundary conditions 

Side boundaries have been assumed to be symmetric and not crossed by any flow. 
At the outer boundary, the flow is generally supersonic and super- Alfvenic; hence, 
information from the outer boundary only propagates downstream, i.e., outside the 
region of interest. In other words, all eight radial characteristic directions at the 
outer boundary are positive. This implies that the boundary conditions at 7 
can be specified arbitrarily. We choose the linear extrapolation method for this 
boundary. 

2.3. Physical boundary conditions 

For the inner boundary, i.e., the physical boundary, the projected characteristics 
method has been used (Hu and Wu, 1984; Wu and Wang, 1987). In the two- 
dimensional problem, two of the six radial characteristic directions are negative 
(t? r is greater than zero but (y r — tv) and (v r — Vf) are smaller than zero, where v s 
and i if are the slow and fast magnetosonic speeds, respectively. v T is the radial flow 
speed); consequently, information from the region of interest propagates upstream 
to the boundary only along these two characteristic directions. In this case, four 
dependent variables at the lower boundary can be specified arbitrarily, and two 
must be calculated from the compatibility relations for every time step. 


3. Pre-Event Coronal Atmosphere 

In order to carry out CME simulations, we first have to define the pre-event corona. 
In the present case, the four pre-event coronal atmospheres have been modeled, via 
a relaxation technique, by prescribing appropriate initial and boundary conditions 
to the set of ideal MUD equations. 

In case 1 we solve the usual set of ideal MHD equations (the reader is referred 
to Wang et al. (1993), for a complete description of the equations), assuming 
constant density (2.25 x 10 8 cm -3 ) and temperature (1.8 x 10 3 * * 6 K) at the solar 
surface, a polytropic Parker-type solar wind solution (with 7 — 1.05) and a dipolar 
potential field configuration with the magnitude 1.67 G at the equator and 3.35 G 
at the pole on the solar surface. The steady-state solutions for density, velocity, 
and magnetic field, shown, respectively, in Figures l(a-c), are identical to those 
of Wang et al. (1993), and are similar to other published models (SH; Guo et al ., 

1991; Mikic and Linker, 1994; and Linker, Van Hoven, and McComos, 1992). 
These solutions adequately represent the streamer, but are unable to reproduce 
the characteristics of the adjacent coronal hole, as easily seen in the behavior of 
the density. Because of this deficiency, Noci et al. (1993) were unable to model 
the observed profile of the La emission vs height in the quiet corona using the 
case 1 model. The present study aims at modeling a pre-event coronal atmosphere 
capable of reproducing simultaneously both the streamer and the hole properties, 
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which may lead to the better understanding of the modeling of both the pre-event 
corona and CMEs. In order to construct such a streamer-hole solution, we adopted 
latitudinal dependent boundary conditions, assuming that the density decreases 
linearly from the equator (1 .43 x 10 8 cm" 3 ) to the pole (1 0 7 cm" 3 ), and temperature 
increases linearly from the equator (1.74 x 10 6 K) to the pole (2.42 x 1() 6 K). 
The magnitude of the magnetic field is 1.36 G at the equator and 2.71 G at the 
pole on the solar surface. We understand that the temperature distribution is not 
realistic. The rationale for this choice lies in the need to provide the polar region 
with the additional energy required to produce a low-density, rapidly accelerating, 
plasma flow. The steady-state solutions obtained by using initial conditions as 
the aforementioned profiles of density and temperature, coupled with the Parker- 
type polytropic solar wind solution (with 7 = 1.05) in a dipolar magnetic field 
configuration, are shown in Figures 2(a-c). These give, respectively, the latitudinal 
variation of the density, of the radial component of the velocity and of the magnetic 
field at several representative heights in the solar corona. Obviously, this solution 
represents both a coronal hole at the pole and a streamer at the equator. 

The choice of 7 being 1.05 for the ideal MHD cases (case 1 and case 2) is based 
on the MHD solar wind flow solution. It assures the solar wind flow at the base is 
subsonic and sub-Alfvenic and reaches supersonic and super-Alfvenic conditions 
at 7 But, the solar wind plasma is a fully-ionized plasma in which case 7 
should be 1.67. Hence, in cases 3 and 4 we adopt 7 = 1.67 and modify the ideal 
MHD energy equation with a source term. To this end, we have added to the energy 
equation a volumetric heating term, similar to the one adopted by SH, given by 

0<-o.i (!i-Nu)/Re) ? ( 1 ) 

where C is an arbitrary constant whose value is dictated by the amount of the 
energy which is being added. Also like SH, we assume that the heat source term is 
independent of latitude and decays exponentially with altitude. 

Figure 3 shows the results we obtain assuming 7 = 1.67, and modifying the 
energy equation as previously described, and using the identical boundary and 
initial conditions as case 1. The results are qualitatively similar to case 1, and do 
not display any of the coronal hole properties because of the constant boundary 
conditions. However, the magnetic field configuration differs from that of case 1, 
as the closed field streamer region shows a sharp cusp reminiscent of observed 
configurations. 

In case 4, we present (Figure 4) the results for latitudinal dependent boundary 
conditions identical to those of case 2; 7 = 1.67, and a modified energy equation 
as in case 3. The latitudinal dependent boundary condition is used here to enable us 
to obtain streamer and hole simultaneously. Like case 3, the closed field streamer 
region shows a sharp cusp. In reference to case 2, we notice that the radial velocity of 
case 4 has increased appreciably because of the added volumetric energy source. 

These four pre-event coronal models will be the initial states wherein CMEs 
will be simulated. Before going on to describe this part of our work, we would like 
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Magnetic Field Lines 



(C) 

Fig. 1 . Latitudinal dependence of density and radial velocity at different heights (a), (b) and magnetic 
field lines (c) for the pre-event corona of case 1. 


to discuss further the pre-event atmospheres. We compute the total energy content 
of the pre-event coronae by integrating four modes (thermal, kinetic, magnetic, and 
gravitational) of energy over the computational domain. These results are given 
in Table I. These results show that, due to the introduction of the additional heat 
source in the energy equation in cases 3 and 4 and the boundary conditions, these 
two pre-event coronae have a higher energy level than cases 1 and 2. 


4. Coronal Mass Ejections in Four Model Atmospheres 

In order to examine the role of pre-event coronae in the modeling of coronal mass 
ejections (CMEs), we applied the same pressure and velocity perturbation at the 
inner boundary of the four model atmospheres, within the closed field region, i.e., 
within the interval 74.25° < 0 < 105.75°. The pressure perturbation is described 
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Magnetic Field Lines 



Fig. 2. Latitudinal dependence of density and radial velocity at different heights (a, b) and magnetic 
field lines (c) for the pre-event corona of case 2* 


TABLE I 


Energy content of the pre-event coronal atmosphere 



Case 1 

Case 2 

Case 3 

Case 4 

Em 

6.00 

4.56 

13.79 

9.55 

Ek 

0.39 

0.66 

2.57 

1.60 

E r 

1.26 

4.05 

8.35 

4.20 

E a 

8.17 

7.45 

13.60 

6.59 

Total 

18.78 

16.72 

38.31 

21.94 


(1() 25 erg km ') 
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Magnetic Field Lines 



7654321 1234567 

(C) 

Fig. 3. Latitudinal dependence of density and radial velocity at different heights (a, b) and magnetic 
field lines (c) for the pre-event corona of case 3. 


by a factor two increase both in the density and temperature; hence, pressure 
increases by a factor of four with respect to its pre-event value. The velocity 
perturbation is set to zero at the edge of the closed field region (i.e., at 74.25°) 
and varies linearly up to a maximum of 200 km s _1 near the equator (i.e., at 87.75 
-92.25°). The velocity also has a linear temporal distribution between / = 0 to 
t = 200 s, defined in time steps of At = 20 s, and the pressure pulse is a step 
function. Both perturbations are maintained throughout the simulation. The total 
energy input resulting from this perturbation amounts to ~ 6 x 10 24 erg km -1 and 
the total mass input is ~ 6.8 x 10 9 g km~ 1 . Taking into account that the latitudinal 
extent over which the perturbation is applied corresponds to seven grid points (or 
3.8 x 10 5 km), the total energy input is 2.4 x 10 30 erg and the total mass input is 
2.6 x 10 15 g. These are comparable to the observed total energy and mass content 
of typical CMEs. 
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Magnetic Field Lines 



Fig. 4. Latitudinal dependence of density and radial velocity at different heights (a, b) and magnetic 
field lines (c) for the pre-event corona of case 4. 


For each of the four pre-event coronae, we give the fractional density enhance- 
ment contours, the magnetic field configuration and the radial velocity at time = 
3000 s and 9000 s, respectively which are summarized in the following. 

Case 1 : dipolar streamer with a polytropic gas, 7 = 1.05 

In this case, the pre-event corona possesses the characteristics of a streamer, with 
no adjacent hole (i.e., no significant density contrast between the pole and equator). 
The density enhancement contours and magnetic field topology for the simulated 
coronal dynamical responses due to the pressure and velocity perturbation in this 
model atmosphere are shown in Figure 5. These results are similar to those given 
by SH, even though they used a perturbation given in terms of an energy input 
within a portion of the closed field region of the streamer together with modified 
energy equation, rather than a boundary perturbation and ideal energy equation as 
used in this study. Both results reproduce a number of observed characteristics of 
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Density Distribution (d-dO)/dO 


Magnetic Field Lines 


t=3000s 




Fig. 5. The evolution of density enhancement (a, c) and magnetic field (b, d) for case 1 at t — 3000 s 
and I = 9000 s after a standard perturbation was applied. The broken lines represent the fast wave 
front and the dotted lines represent the density depletion. 


loop transient, including: (1) the high-density region, in the shape of Megs’, and 
(2) the rarefaction region (low density) between the two legs. The propagation 
speed measured from the apex of the loop is ~ 250 km s” 1 which corresponds 
to a typical slow CME speed. This propagation speed is found to depend on the 
strength of the perturbation at the lower boundary. 

Case 2: dipolar streamer and coronal hole with a poly tropic gas , 7 =1.05 
The simulated density enhancement contours and magnetic fields for this pre-event 
corona are shown in Figure 6. These results are basically similar to case 1 except 
that the loop system exhibits a larger latitudinal extension than case 1. This is due 
to the accompanying low-density coronal hole which has a larger fast-mode wave 
speed. From Figure 7 the initial ft distribution in case 2 shows that when the density 
decreases from the equator to the pole the ft value varies in latitude from O[10] to 
O[10 -1 ]. In case 1, the density is almost uniform in latitude, thus the ft value along 
latitude lines is also almost constant. By comparing the fast mode speeds in these 
two cases, the fast speed in case 2 along a latitude line is larger than that in case 1. 
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Density Distribution (d-dO)/dO 


Magnetic Field Lines 


t=3000s 




(a) 


(b) 


t=9000s 



Fig. 6. The evolution of density enhancement (a, c) and magnetic field (b, d) for case 2 at t = 3000 s 
and t = 0000 s after a standard perturbation was applied. The broken lines represent the fast wave 
front and the dotted lines represent the density depletion. 


Hence, the perturbation propagates more rapidly in latitude in case 2 than in case 1. 
The radial propagation speed for this case is ~ 200 km s " 1 which is lower than in 
case 1 , possibly because in this case loops show some latitudinal expansion. 

Case 3: dipolar streamer in a heated atmosphere with 7 — 1.67 
In this case, the density enhancement contours and magnetic field topology are com- 
pletely different from those obtained in the two previous simulations, as shown 
in Figure 8 . These results show two-stages of development. In the early stage 
(t < 3000 s), the simulation shows a low-lying density and magnetic loop config- 
uration (Figures 8 (a) and 8 (b)), while the magnetic field topology shows both loop 
shape and a magnetic island at t = 6000 s (Figures 8 (c) and 8 (d)). In the later stage 
(t > 6000 s), the magnetic island forms a distinct plasmoid shape which moves 
outward at a speed of 275 km s -1 with a loop shaped density enhancement and 
magnetic field lines extending outside the computational domain in Figures 8 (e) 
and 8 (f). In this case the pre-event corona is in the highest energy level in compar- 
ison with the other three cases. From Figure 7(c) the 3 values at the equator near 
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Pre-event Plasma Beta Distribution 


Case 1 Case 2 



Case 3 



Fig. 7. The two-dimensional distributions of plas 
respectively. 


Case 4 



the solar surface within a small region are two order of magnitude larger than in the 
polar region. Thus, when velocity and pressure pulses were introduced the density 
in this region easily increases because the plasma is easily compressed. The field 
lines are pushed together to induce a pinch effect and to change the field topology 
into that believed to exist in magnetic clouds (Suess, 1988). 

Case 4: dipolar streamer and coronal hole in a heated atmosphere with 7 = 1.67 
The results for this case (shown in Figure 9) show that there are two bright legs 
(density enhancements) extending to the high corona, forming a spike-like struc- 
ture with an overlying loop. According to the description of CME morphologies 
(Burkepile and St. Cyr, 1993), case 4 more like the mound: ‘the tops of mounds 
often have a well-defined, curved appearance similar to the frontal loop, but there is 
no obvious decrease in brightness behind the the leading edge’. Unlike case 3, the 
magnetic field lines do not form a plasmoid. In this case the energy level initially 
is almost a factor two smaller than in case 3, but much larger than cases 1 and 2. 
Further, as shown in Figure 7(d), the fi distribution is similar to case 2 in which the 
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Density Distribution (d-dOydO 


Magnetic Field Line* 


txJOOOe 


(*> M 


U6000s 


(c) (d) 


t* 9000s 



and / = 9000 s after a standard perturbation was applied. The broken lines represent the fast-wave 
front and the dotted lines represent the density depiction. It should be noted that the broken line does 
not appear at / = 9000 s, because the fast waves have already propagated out of the computational 
domain. 






value is smaller at the pole and larger at the equator, because the density is smaller 
at the pole and larger at the equator. But Figure 7 shows that the ft values in case 4 
(Figure 7(d)) are much smaller than other three cases. This will lead to much easier 
propagation of the disturbances in the radial direction. Also the small ft means that 
the magnetic field strength is strong, thus, the field lines are not easy to compress 
by latitudinal flow at the equatorial region to cause the formation of a plasmoid as 
in case 3. 
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Density Distribution (d-dO)/dO 


Magnetic Field Lines 




Fig. 9. The evolution of magnetic field and density enhancement for case 4 at t — 3000 s and 
t = 9000 s after a standard perturbation was applied. The broken lines represent the fast-wave front 
and the dotted lines represent the density depletion. It should be noted that the broken line does 
not appear at f = 9000 s, because the fast waves have already propagated out of the computational 
dom ain. 


5. Discussions and Concluding Remarks 

We have presented results from a numerical simulation of CMEs in four different 
pre-event coronae in which two cases possess both characteristics of streamer and 
hole, and have shown how pre-event coronal models affect CMEs simulations. In 
this section we discuss the new features of our work and compare it with previous 
numerical simulations of CMEs. 

From these simulations, we obtain three distinct types of CMEs: (i) classical 
loop-like transients, (ii) plasmoid enveloped by a loop, and (iii) spike-like legs in 
magnetic field configuration and mound shape of density enhancement. All these 
results are obtained using a standard perturbation with four different pre-event 
coronae. Now the question which needs to be answered is why? Consider the 
energy content of the pre-event coronae, as given in Table I. This shows that case 3 
has the highest energy content of all models. The reason for the variation of energy 
content is the different boundary conditions and additional energy being added in 
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Plasma Beta Distribution at 9000s 


Case 1 Case 2 



Case 3 Case 4 



Fig. 10. The two-dimensional distributions of plasma p for evolution states at 9000 s of case 1, 2, 
3, and 4, respectively. 


cases 3 and 4. Although we are unable, at this time, to give threshold values for the 
occurrence of a plasmoid, we surmise that a plasmoid-type CME occurs only when 
the pre-event corona has a sufficiently high-energy content. The middle range of 
energy content is responsible for mound-type CMEs, and the lowest range gives the 
classical loop-like CMEs. This conjecture coincides with observations reported by 
Burkepile and St. Cyr (1993), who ranked the morphologies of CMEs according to 
their relative frequency of appearance, showing that a loop-like transient is ranked 
first, followed by mound shape (in the fifth place) and blobs last (seventh place). 

Steinolfson (1988) suggested that loop-like CMEs would originate in a pre- 
event atmosphere, where shocks would not easily form. Simulations based on the 
present pre-event coronal models allow us to confirm this claim. In the four cases 
here the highest density regions are behind the fast wave propagation lines. As 
shown, for instance, in Figure 9, the MHD fast-mode speed increases rapidly from 
the equator to the pole, along the meridional direction, as well as from the solar 
surface up to the open field regions, along the radial direction, thus preventing the 
formation of fast shocks. 




